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ABSTRACT 

A  single-stage  turbine  is  designed  by  utilizing  the  circula¬ 
tion  control  method  [1],  This  method  has  broken  the  conventional 
rules  of  design.  A  model  experiment  is  conducted  in  order  to 
verify  the  reliability  of  the  method.  General  characteristics  are 
measured  over  a  wide  range.  The  measurements  of  the  efficiency 
distribution  along  the  blade  span  and  of  the  inlet' and  exit  velo¬ 
city  triangles  of  the  rotor  blades  are  made  around  the  design 
point.  The  main  test  results  are  given  in  this  paper  and  a  simple 
analysis  is  made.  The  results  Indicate  that  the  main  characteristics 
istics  of  the  circulation  control  design  method  are  verified  and 
the  design  requirements  are  satisfied.  The  method  is  practicable 
basically . 

I.  INTRODUCTION 


The  design  of  aircraft  turbines  has  been  conducted  according 
to  traditional  equi-circulation,  equal  to  aj  or  the  mixed  rules  derived 
from  them.  The  variation  of  the  degree  of  reaction  along  the 
blade  span  is  relatively  large.  Numerous  experimental  investiga¬ 
tions  have  shown  that  the  flow  conditions  at  the  root  and  the  tip 
are  not  desirable.  The  second  order  loss  and  the  leakage  loss  at 
the  tip  of  the  blade  is  relatively  large.  With  the  advancement 

1 - 
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made  in  the  study  of  the  internal  flow  of  turbines  and  in  the 
application  of  computers,  it  is  possible  now  for  the  designers  to 
bypass  the  traditional  design  laws.  Utilising  the  complete  radial 
equilibrium  equation  to  control  the  distribution  of  the  circula¬ 
tion  V^R  along  the  blade  span,  the  degree  of  reaction  at  the  root 
is  raised  while  the  degree  of  reaction  at  the  tip  is  lowered. 

The  flow  conditions  at  both  the  root  and  the  tip  are  improved  and 
so  is  the  performance  of  the  turbine.  In  1972  we  established  a 
research  group  to  conduct  work  in  this  area.  A  computer  program 
was  developed  in  coordination  with  the  computing  laboratory  at  the 
Academy  of  Sciences  and  calculations  were  made.  In  May  1975,  a 
single-stage  turbine  was  designed.  This  is  the  first  trial  violat¬ 
ing  the  conventional  design  method.  Reference  [1]  describes  the 
method  of  circulation  control  design  and  its  main  characteristics. 
The  main  results  of  the  design  are  also  presented.  In  order  to 
verify  the  reliability  of  this  new  design  method  and  to  draw  a  con¬ 
clusion,  a  test  model  was  designed  and  a  first  stage  testing  is 
carried  out  on  the  model  turbine  tester  in  Shenyang  from  April  to 
July,  1978.  In  this  paper,  the  principal  results  of  the  test  are 
presented  and  the  fundamental  analysis  is  made. 

1.  Test  model 


The  prototype  is  a  single-stage  turbine.  When  designing,  the 
complete  radial  equilibrium  equation  Is  employed  at  the  gaps  in  the 
axial  direction  to  control  the  distribution  of  the  circulation  VQR 
along  the  blade  span,  intending  to  improve  the  flow  conditions  at 
the  root  and  the  tip.  The  work  along  the  blade  span  is  linearly 
distributed  on  the  whole.  The  variation  of  the  work  is  7-56%. 
Uniform  velocity  coefficient  distribution  is  takeh  along  the  blade 
span,  that  is,  for  the  guiding  blade,  $  =  0.96  and  for  the  rotor 
blade,  =  0.95.  The  design  load  is  relatively  high.  The  main 

parameters  are:  effective  work  Lm  =  21,770  kg-m/kg,  loading  coeff- 

„ 1  * 
icient  at  the  midspan  is  P  a  2.085,  stagnation  expansion  ratio  % 

ft  ^ 

*  2.39,  design  efficiency  nT  =  0.879.  The  Mach  number  at  the  exit 
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of  the  guiding  blades  ’■  r.  relatively  high  and  the  flow  Is  found  to 
be  supersonic  near  the  root  region.  At  the  root  Mc^  =  1.11^.  The 
scale  of  the  test  rr.odel  is  m  =  0.75*  The  flow  channels,  the  main 
measuring  positions  and  parameters  are  shown  in  Figure  1. 

The  relative  measurement  error  of  the  stagnation  efficiency, 
the  relative  net  flow  rate,  and  the  net  work  at  the  large  expansion 
ratio  condition  is  around  ±  1.02. 

The  rotor  blades  are  capped  and  at  the  cold  condition,  the 
radial  gap  at  the  blade  tip  is  0.5-0. 8  mm,  which  is  0.8-1.27%  of 
the  blade  span. 


Figure  1.  The  flow  channel  and  the  measuring  positions 
2.  Test  results  and  analysis 

(a)  Performance  at  the  design  point 

The  tested  values  and  the  designed  values  of  the  main  perfor¬ 
mance  parameters  are  compared  below: 

(b)  Stagnation  efficiency 

Based  on  the  average  value  of  the  multiple  points  total  press¬ 
ure  at  section  2-2  far  downstream  and  the  torque  of  the  hydraulic 

* 

work  detector,  the  stagnation  efficiency  nTQ  2  of  the  single-stage 
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Stagnation  efficiency  I  net 


measured 
value  at 
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stream 
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nT0-2 
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measured 
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value 


designed 
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net  flow 

of  the  guid- 
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ing  blade  , 
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4 

0.99 

229 

1.0 

232.8 

12 

-1.652 

turbine  is  calculated.  The  efficiency  at  the  design  point  meets 
the  original  design  requirement. 


It  should  be  pointed  out  that  the  efficiency  calculated  by  the 
single-point  total  pressure  at  near  downstream  is -1. 5-^2  higher 
than  that  calculated  by  the  multiple  points  total  pressure  at  far 
downstream.  The  deviation  of  the  expansion  ratio  is  relatively 
larger,  which  is  2.8%  higher  than  that  at  the  design  point.  After 
analyzing,  it  is  more  desirable  to  take  the  far  downstream  results. 
The  static  probe  behind  the  rotor  blade  detects  a  pulsive  flow  with 
a  frequency  of  about  10,000  Hz.  When  the  tangential  velocity  is 
larger  than  the  absolute  velocity  of  the  flow,  the  pressure  detected 
is  higher  than  the  true  value.  The  existence  of  a  vortex  flow 
behind  the  rotor  blade  leads  to  a  loss  in  the  total  pressure.  The 
loss  cannot  be  measured  accurately  at  the  near  downstream  since  the 
probe  suffers  the  "choking  effect"  if  it  is  placed  close  to  the 
trailing  edge.  Hence,  In  order  to  make  accurate  measurements,  there 
should  be  a  straight  section  behind  the  rotor  blade.  The  fixed 
probe  should  be  placed  at  the  nearly  uniform  region  of  the  flow.  For 
the  transonic  turbine,  this  region  Is  two  or  three  cascade  widths 
behind  the  cascade  [2],  At  the  point  122.8  mm  away  from  the  trail¬ 
ing  edge  of  the  guiding  blade,  there  is  still  a  wave  motion  of  37  nun 

Hg  column  (Figure  2).  This  wave  is  induced  by  the  wake  of  the  guid- 
,1 

ing  blade  passing  through  the  rotor  blade.  This  phenomenon  indicates 


that  the  wake  of  the  transonic  cascade  is  relatively  long  and 
that  it  is  more  desirable  to  calculate  the  total  pressure  by  the 
static  pressure  at  the  exit  and  other  measured  parameters  [3]. 


(i  :'i(,  k  in 
^  •'/itfcK'AiHO.ftM  Win  'i  Ht-A  1 M  *> 


n(No.  of  revolutions  of  slide  plate,  each  revolution  corresponds 
to  a  movement  of  1.5  mm). 

Fig.  2.  The  variations  of  the  total  pressure  and  the  static 
pressure  along  the  cascade  width  far  downstream  of  the  rotor  blade, 
(conditions  are  maintained  with  an  error  of '4  mm  Hg). 

(c)  The  relative  net  flow  rate 

It  Is  evident  from  the  tested  characteristics  that  the  throat 
of  the  guiding  blade  at  the  design  point  is  the  critical  cross- 
sectional  area.  The  guiding  blade  throat  area  of  the  test  model 
is  1.65$  smaller  than  the  designed  value,  causing  the  net  flow 
rate  to  be  slightly  less  than  the  designed  value. 

II.  INTERNAL  FLOW  PATTERN  OF  THE  TURBINE  NEAR  THE  DESIGN  POINT 

In  order  to  study  the  Internal  flow  pattern  in  the  turbine  of 

circulation  control  design,  measurements  along  the  blade  span  are 

# 

made  near  the  design  point  (n  =  1.0,  tt^,  =  2.3).  Measurements  are 
made  with  the  combination  probe  for  the  total  pressure,  the  direc¬ 
tion,  the  static  pressure  and  total  temperature  at  the  nine  radial 
positions  along  the  ^5°  line  19.5  mm  from  the  guiding  blade  and 
76  mm  from  the  rotor  blade.  It  is  also  moved  along  the  cascade 
width  behind  the  guiding  blade.  The  results  are  tabulated  through 
a  wire-expansion  type  transducer  and  recorded  continuously  by  an 
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automatic  six  point  recorder.  The  mean  value  is  calculated  accord¬ 
ing  to  the  area  integral.  For  the  region  where  m  >  1,  the  total 
pressure  and  the  linearly  distributed  values  of  the  static  press¬ 
ure  on  the  inner  and  outer  walls  detected  by  the  probe  are  cor¬ 
rected  according  to  the  normal  shock  relations. 

1 .  Flow  parameters  at  the  exit  of  the  guiding  blade 

One  essential  characteristic  of  the  circulation  control  design 
is  the  "reverse  twist"  of  the  exit  angle  of  the  guiding  blade,  that 
is,  the  angle  is  greater  at  the  root  and  smaller  at  the  tip.  The 
measured  value  along  the  blade  span  is  relatively  smooth  and  is 
almost  parallel  to  arcsin(a/t)  of  the  guiding  blade  except  at  both 
ends.  If  a  translation  is  made,  it  is  observed  that  they  are  iden¬ 
tical  for  most  of  the  region.  The  deviation  is  due  to  systematic 
error.  The  measured  values  and  the  trend  of  the  exit  angle  do 

not  agree.  Calculation  utilizing  the  measured  values  indicated 
that  the  ^  at  the  root  is  larger  than  the  designed  value.  These 
results  are  shown  in  Figure  3* 

A  linear  distribution  of  static  pressure  between  the  inner  and 
outer  walls  is  taken  at  the  exit  of  the  guiding  blade.  The  exit 
total  pressure  and  the  exit  velocity  coefficient  along  the  blade 
span  are  almost  constant  for  most  part  of  the  region  except  at  the 
two  ends,  where  they  are  slightly  smaller  (Figure  3,  Figure  5). 

This  indicates  that  the  second  order  loss  of  the  guiding  blade  is 
not  large  at  all  and  the  region  affected  is  relatively  small.  The 
mean  value  of  the  velocity  coefficient  $0_1  =  0.9815  is  higher  than 
the  designed  value  0.96.  This  is  because  the  measurements  made 
do  not  detect  the  loss  completely.  The  flow  afterwards  will  still 
be  undergoing  a  mixing  process. 

Another  characteristic  of  the  circulation  control  design  is 
the  gradual  increase  of  the  slope  and  the  curvature  of  the  stream¬ 
lines.  The  streamlines  are  squeezed  against  the  root.  This  is 
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calculated  from  the 


proved  by  the  relatively  large  value  of  M  ^ 
measured  results.  This  also  justifies  the  necessity  of  consider¬ 
ing  the  slope  and  curvature  of  the  streamlines  in  the  circulation 
control  design  (Figure  3). 


isured  value 
— x — arc  sin  ( a/t ) 

~*~alk 

— x — designed  value 

*  2 
p  (kg/cm  ) 


Figure  3.  The  distribution  of  the  exit  parameters  of  the 

guiding  blade  along  the  blade  span. 

2 .  Flow  parameters  at  the  exit  of  the  rotor  blade 

Measurements  along  the  blade  span  are  made  at  the  far  down¬ 
stream  section  2-2.  The  measurement  for  the  absolute  flow  angle  at 
the  exit  is  repeated  twice  and  is  found  to  be  consistent.  There  is 
a  slight  deviation  in  the  lower  half  region  but  its  value  is  closer 
to  the  designed  value.  Within  the  region  j  blade  height  from  the 
tip,  the  deviation  is  larger  (up  to  20°)  due  to  the  leakage  loss 
effect.  Comparison  between  the  relative  flow  angle  and  the  arcsin 
(a/t)  of  the  rotor  blade  shows  a  similar  situation.  The  deviation 
at  the  tip  is  approximately  5°-  (Figure  4). 

The  exit  total  pressure  is  relatively  low  at  the  mid  span  and 
is  relatively  high  at  the  30%  and  80%  blade  height.  It  falls  sharnly 
at  both  ends,  producing  a  saddle  shape.  The  exit  static  pressure  is 
almost  linearly  distributed  along  the  blade  height  (Figure  4).  The 
off-design  does  not  cause  any  large  non-uniformity  of  the  exit 
parameters . 
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Figure  A.  Distribution  of  the  measured  parameters  far 
downstream  of  the  rotor  blade  along  the  blade  span 


With  the  inlet  and  exit  parameters  obtained,  the  rotor  blade 
velocity  coefficient  ^1_2  is  calculated,  which  reaches  0.97  in  the 
mainstream  region.  It  falls  sharply  at  both  ends  .since  the  inten¬ 
sity  and  the  affected  region  are  both  increased  with  respect  to  that 
of  the  guiding  blade  due  to  the  secondary  flow  effect  (Figure  5). 

Its  mean  value  is  approximately  0.91>  which  is  much  lower  than  the 
designed  value  0.95.  This  velocity  coefficient  is  obtained  from 
the  measured  parameters  far  downstream.  The  mixing  loss  at  76  mm 
behind  the  rotor  blade  Is  considered,  and  hence  the  measurement  of 
the  loss  is  more  complete.  Since  the  wake  from  the  guiding  blade 
is  comparatively  long,  it  is  still  under  a  mixing  process  when  it 
reaches  the  rotor  blade  and  even  behind  the  rotor  blade.  This 
causes  an  undesirable  operating  environment  for  the  rotor  blade. 
Hence,  part  of  the  loss  detected  is  induced  by  the  guiding  blade. 
During  the  design  process,  this  kind  of  practical  situation  should 
be  dealt  with. 


3.  Stagnation  efficiency 

Based  on  the  total  temperature  and  the  total  pressure  at  the 

« 

Inlet  and  exit,  the  stagnation  efficiency  ru  is  computed.  At  the 

ft  ^ 

75 %  blade  height,  the  efficiency  reaches  a  maximum  of  0.95.  Within 
the  region  30%  blade  height  from  the  root,  the  decrease  of  the 
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Figure  5.  The  distribution  of  Figure  6.  Variation  of  the 
velocity  coefficient  of  the  guid-  stagnation  efficiency  along  the 
ing  blade  ^  and  the  velocity  blade  span, 

coefficient  of  the  ro: or  blade 

ip  along  the  blade  span.  — o---experimental  value 

(n  =  1.0,  ttJ,  =  2.301) 

- designed  value 

(n  =  1.0,  7T*  =  2.39) 

efficiency  is  large.  It  also  falls  within  the  region  20%  blade 
height  from  the  tip  (Figure  6).  Employing  the  variation  design, 
the  function  of  the  tangential  velocity  is  enhanced  and  a  favorable 
effect  on  the  efficiency  is  produced. 


From  the  results  and  discussion  presented  above,  the  main 
characteristics  of  the  control  circulation  design  are  verified  and 
the  original  design  requirements  are  satisfied  on  the  whole.  The 
secondary  flow  behind  the  guiding  blade  is  not  so  destructive. 
However,  the  effect  due  to  the  secondary  flow  behind  the  rotor  blade 
and  the  leakage  loss  at  the  tip  is  quite  significant.  This  causes 
a  relatively  large  deviation  of  the  flow  parameters  from  the 
designed  values  at  both  ends,  especially  the  deviation  of  the  abso¬ 
lute  exit  flow  angle  at  the  blade  tip  (20°).  Further  consideration 
on  the  variable  loss  design  should  be  made.  This  does  not  only 
require  the  second  order  leakage  loss  model  at  the  tip,  but  also  the 
model  concerning  the  effect  of  the  exit  flow  angle.  Further  studies 
on  all  these  models  are  required. 
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Figure  7.  The  stagnation  efficiency  vs  and  n 


III.  TEST  FOR  GENERAL  CHARACTERISTICS 


For  the  stagnation  expansion  ratio  ti^  =  1.3-3.08,  and  the 
converted  speed  of  rotation  n  =  0.5-1 *1>  the  general  characteris¬ 
tics  arc  tabulated  (Figure  7>  Figure  8,  Figure  9).  How  these 
characteristic  lines  vary  from  the  general  rules  is  not  important 
and  will  not  be  discussed  in  detail  here.  However,  the  following 
problems  should  be  discussed: 


Figure  8.  _The  relative  net 
flow  rate  G  vs  tt* 


1. 


The  internal  operating  condition  of  the  turbine  under 


different  states 
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The  internal  operating  condition  of  the  turbine  at  the  off- 

design  condition  can  be  observed  clearly  from  the  variation  of  the 

guiding  blade  expansion  ratio  p*/p^,  the  rotor  blade  expansion 

ratio  p^/pi,  and  the  stagnation  expansion  ratio  pVpi,  vs  the  state 

(Figure  10).  The  corresponding  curve  of  the  degree  of  reaction  is 

also  shown  (Figure  11).  Degree  of  reaction  p?  T  ill 

„  i  -  ip: /pi )-«*-*v* 

With  the  increase  of  tt*,  the  guiding  blade  expansion  increases 

linearly.  Over-expansion  appears  at  the  root,  causing  a  pressure 
diffusion  condition  (p^/pi,  <  1,  PT  <  0)  at  the  root  of  the  rotor 
blade.  The  root  of  the  guiding  blade  reaches  the  critical  condition 
quickly  and  the  inclined  plane  exit  expansion  appears  right  after¬ 
wards.  As  the  blade  reaches  the  critical  condition,  the  pressure 
perturbation  behind  it  cannot  be  transmitted  upstream.  This  prevents 
the  inclined  plane  exit  expansion  at  the  root  of  the  guiding  blade 
from  reaching  the  limit.  Since  the  tip  of  the  guiding  blade  does 
not  reach  the  critical  point,  the  pressure  there  will  no  longer 
increase.  However,  the  inclined  plane  exit  can  still  be  expanding 
until  its  limit  expansion  ratio  is  reached.  Hence  the  total  expan¬ 
sion  ratio  increases  continuously  and  the  corresponding  degree  of 
reaction  increases  linearly.  From  the  trend  of  the  expansion  ratio 
curve,  it  can  also  be  seen  that  the  circulation  control  design  makes 
the  flow  in  the  radial  direction  to  become  more  gradual.  This  does 
not  only  favor  the  design  states,  but  also  improve  the  non-designed 
states,  especially  the  performance  of  the  total  expansion  ratio. 

The  variations  of  the  expansion  ratio  and  the  degree  of  reaction 
with  the  speed  of  rotation  are  more  or  less  similar. 

2 .  The  oscillation  of  the  stagnation  efficiency  at  the 
transonic  rerion 


Oscillation  of  the  efficiency  curve  appears  near  the  design 
point.  A  relatively  greater  number  of  measuring  points  is  taken 
here  and  measurements  are  repeat  a  ''nee.  These  readings  are  found 
to  be  consistent,  indicating  that  it  is  not  simply  an  accidental 
phenomenon.  The  reasons  are  analysed  below:  When  designing  the 
expansion  ratio  of  the  turbine,  the  operating  condition  of  the 


guiding  blade  la  unchanged  basically,  and  so  Is  the  corresponding 
angle  of  attack  of  the  rotor  blade.  The  exit  Mach  number  of  the 
rotor  blade  is  in  the  range  0.8-1. 0.  Hence,  shockwaves  are  liable 
to  occur  at  the  throat  and  the  trailing  edge.  The  shockwaves  from 
the  trailing  edge  extend  to  the  neighboring  blade  and  cause  separa¬ 
tion  of  the  boundary  layer.  When  tt*  increases  the  shockwaves  lean 
back  and  the  separation  point  moves  back.  As  the  blade  profile 
becomes  slender  and  the  second  order  loss  decreases,  n*  further 
increases.  The  intensity  of  the  shockwaves  increases  and  hence  the 
loss  increases,  producing  an  oscillating  phenomenon.  This  pheno¬ 
menon  is  common  for  plane  cascade  and  ring  cascade. 


Figure  10.  The  variation  of  the  expansion  ratio  vs  state 


1? 


,  / 


designed  value: 


Figure  11.  Variation  curve  of  the  thermodynamic  degree  of 
reaction  vs  state 


Transonic  design  is  an  effective  way  to  increase  the 
stage  load 


From  the  curves  of  the  stagnation  efficiency -and  the  net  work, 
it  can  be  observed  that  the  decrease  of  the  efficiency  curve  is 
not  very  sharp  at  relatively  high  it*.  But  there  is  still  a  large 
growth  of  the  stage  load.  The  related  parameters  are  listed  in 
the  table  below: 


”T 

2.1 

2.6 

2.7 

2.8 

2.85 

2.9 

3.0 

Vt 

0.881 

0.883 

0.876 

0.861 

0.86 

0.852 

0.81 

It 

20 

21.5 

22.1 

22.5 

22.6 

22.7 

22.8 

^’itbiUVo 

0 

—  0.113 

-0.905 

-2.26 

-2.71 

-3.5 

-1.86 

Alt/Letfl 

0 

7.5 

10.5 

12.5 

13 

13.5 

11 

The  variation  of  the  efficiency  for  it#  in  the  range  2.  *1-2. 7 
is  relatively  small,  while  the  net  work  increases  almost  linearly. 
The  efficiency  curve  decreases  linearly  afterwards’.  The  increase 
of  the  net  work  gradually  decelerates  since  expansion  ratio  tends 
to  reach  its  limit.  As  the  expansion  ratio  increases  to  2.7,  the 
efficiency  falls  only  0.905%  while  the  net  work  increases  by  10!!. 

At  this  point,  the  relative  exit  Mach  number  at  the  mid  span  is 
estimated  as  1.05.  The  main  factor  affecting  the  efficiency  in  this 
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region  is  the  exit  Mach  number.  Hence  the  transonic  design  is  a 
prospective  way  of  increasing  the  stage  load  without  greatly  reduc¬ 
ing  the  efficiency.  Process  in  this  area  can  be  expected  through 
certain  solid  theories  and  experimental  work. 

IV  .  SU i-li” ARY 

The  following  can  be  concluded  from  the  above  analysis: 

1.  The  turbine  is  designed  by  the  circulation  control  method. 
The  principal  performance  parameters  at  the  design  point  satisfy 
the  design  requirement  basically.  The  measurements  along  the  blade 
span  near  the  design  point  show  that  the  principal  characteristics 
of  the  circulation  control  design  method  are  verified.  Hence  the 
method  is  applicable  basically. 

2.  The  secondary  flow  behind  the  guiding  blade  is  not  too 
essential  while  the  intensity  of  the  secondary  loss  behind  the  rotor 
blade  and  the  region  affected  increases.  Due  to  the  leakage  at  the 
blade  tip,  the  exit  flow  angle  at  the  tip  deviates  greatly  from  the 
design  value.  Hence,  further  consideration  of  variation  loss  design 
is  necessary. 

3.  The  expansion  ratio  curve  tabulated  by  means  of  the  sta¬ 
tic  pressure  at  the  wall  is  very  favorable  to  the  understanding  of 
the  flow  process  inside  the  turbine.  The  analysis  shows  that  the 
circulation  control  design  does  not  only  raise  the  performance  of 
the  design  state,  but  also  improves  the  performance  of  the  non-design 
states,  especially  the  total  expansion  ratio. 

*1 ,  As  the  expansion  ratio  increases  from  2 .  *1  at  the  design 
point  to  2.7,  the  efficiency  decreases  by  approximately  155,  while 
the  effective  work  increases  by  approximately  1055  .  This  indicates 
that  the  trapsonic  stage  design  for  relatively  high  exit  Mach  num¬ 
ber  is  a  very  prospective  way  to  increase  the  stage  load  without 
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reducin':  the  efficiency  significantly. 

5.  The  efficiency  measured  at  the  near  downstream  recion  of 
the  rotor  blade  at  the  desirn  point  is  2.8*  higher  than  that  mea¬ 
sured  at  the  far  downstream  recion.  Hence,  the  section  where 
measurements  are  made  should  be  placed  at  the  far  downstream  of 
the  rotor  blade  exit,  where  the  flow  is  more  uniform. 
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THE  EXPERIMENT  RESEARCH  ON  THE  TURBINE  DESIGN  OF 
CIRCULATION  CONTROL  METHOD 

Group  for  Turbine  Design  in  Designing  Department  of  Zhuzhou  Aircraft 
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Abstract 

A  single-stage  turbine  is  designed  by  using  a  circulation  control  method.  It  has 
broken  the  conventional  law  of  design.  A  model  experiment  was  conducted  in  order  to 
demonstrate  the  reliability  of  this  method.  A  general  characteristic  was  measured  over 
a  wid  range.  The  efficiency  distribution  along  blade  span  and  the  inlet  and  outlet 
velocity  triangle  measurement  of  rotor  blades  were  made  around  the  designing  point. 
The  main  test  results  are  given  and  the  simple  analysis  is  applied  to  it.  The  results  in¬ 
dicated  that  the  main  characteristics  of  . . .  control  design  got  demonstra¬ 

tion  and  the  design  retirement  has  been  reached  on  the  whole.  This  method  is  practi¬ 
cable.  ’ 
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THE  CALCULATION  CF  TRANSONIC  TURBINES 
AT  OFF-DESIGN  CONDITION 


Ge  Man-Chu 


Institute  of  Engineering  Thermophysics ,  Academia  Sinica 


Abstract 


A  group  of  simultaneous  equations,  including 
S0  stream  surface  equations,  simplified  stream 

surface  equations  (for  solving  the  triangular  region 
at  the  tail  of  the  cascade),  flow  loss  correlation 
and  deviation  angle  correlation  are  obtained  for 
the  calculation  of  the  direct  problem  of  turbo- 
machinery.  By  the  solution  of  these  equations,  the 
design  point  as  well  as  the  off-design  performances 
can  be  calculated. 

The  essential  difference  in  treating  subcriti- 
cal,  critical,  and  supercritical  flow  condition  is 
considered.  The  influence  of  losses  is  taken  into 
account  in  the  calculation  of  choking  conditions.  A 
unified  approach  is  adopted  to  solve  these  equations. 
Losses,  deviation  angles,  and  other  flow  parameters 
are  obtained  in  the  same  iterative  process,  so  that 
false  choking  conditions  will  not  occur  and  conver¬ 
gence  will  be  rapid .  The  computer  program  is  then 
completed  and  may  be  used  in  the  calculation  of  the 
off-design  performance  of  transonic  turbines. 

Satisfactory  results  are  obtained  when  applying 
this  program  to  calculate  the  exit  performance  of 
the  turbine  stage.  The  reduced  flow  rate,  the  effi¬ 
ciency  and  the  expansion  ratio  are  fairly  consistent 
with  the  experimental  results.  The  difference  between 
the  calculated  results  and  the  experimental  data  is 
within  1%. 


I.  INTRODUCTION 


In  this  paper,  the  stream  surface  equations  of  the  3-D 
flow  theory  of  turbomachinery,  the  simplified  surface  of 
revolution  equations  (for  the  triangular  region  at  the  tail  of 
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the  cancati"),  the  Ion::  correlation  equation  and  the  equation 
for  the  flow  ancle  at  the  exit  of  the  cascade  are  applied 
simultaneously  to  solve  the  direct  problem  of  turbomachinery . 

The  calculations  of  the  losses  and  the  exit  flow  ancle 
under  different  operating  conditions  are  made  with  this  set 
of  simultaneous  equations.  Hence  not  only  the  flow  parameters 
and  the  performance  at  the  design  point  can  be  calculated, 
but  the  flow  parameters  and  performance  at  off-design  conditions 
may  also  be  calculated. 

The  different  characteristics  of  subcritical,  critical, 
and  supercritical  flow  are  taken  into  account  during  the  cal¬ 
culation  process.  Corresponding  consideration  of  the  losses 
and  exit  flow  angle  is  also  made.  At  the  same  time,  the  influ¬ 
ence  of  the  losses  on  the  calculation  of  the  choking  condition 
is  also  considered.  Hence  a  more  reasonable  result  for  the 
choking  flow  rate  can  be  obtained  readily.  For  the  calculation 
of  the  losses  and  the  exit  flow  angle,  a  unified  approach  to 
the  set  of  simultaneous  equations  is  applied  to  obtain  the 
losses,  flow  angle,  and  other  flow  parameters  at  the  same  time. 
This  approach  prevents  the  occurrence  of  the  false  choking 
condition  which  is  due  to  the  iteration  process  for  the  losses 
and  the  exit  angle.  These  parameters  are  found  to  be  in  good 
synchronism  during  the  calculation.  This  is  favorable  to  the 
convergence  of  the  solution  of  the  set  of  simultaneous  equa¬ 
tions,  and  the  convergence  speed  is  greatly  accelerated. 

A  computer  program  for  calculating  the  transonic  turbine 
at  off-design  conditions  is  developed  and  is  executed  on  the 
TQ-16  machine.  For  a  single  stage  turbine,  the  computation 
time  for  each  point  is  around  20  seconds.  Calculations  are 
made  for  certain  domestic  as  well  as  foreign-designed  turbines. 
For  Instance,  off-design  calculation  is  conducted  for  the  250 
kw  gas  engine  turbine  designed  and  tested  by  the  Institute  of 
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Mechanics,  Academia  fdnic.a,  and  fairly  satisfactory  results 
are  obtained.  For  the  entire  test  range,  the  calculated  values 
of  the  reduced  flow  rate,  efficiency  and  expansion  ratio  are 
found  to  be  consistent  with  the  experimental  data.  The  error 
is  only  about  1G.  A  trial  calculation  is  also  conducted  for 
the  turbine  of  the  JT3D-3B  engine,  and  it  is  found  that  the 
deviation  of  the  calculated  values  of  the  three  parameters  from 
the  given  design-point  data  is  also  around  1% . 

II.  CALCULATION  OF  TRANSONIC  OFF-DESIGN  CONDITION 

To  meet  the  development  of  both  defense  and  civil  indus¬ 
tries,  the  performance  of  the  turbine  at  off-design  conditions 
is  gradually  considered  as  a  design  parameter  also.  The  two- 
element  method  is  gradually  employed  to  calculate  the  perform¬ 
ance  of  the  turbine  at  off-design  conditions.  But  there  is  a 
deficiency  in  the  calculation  of  losses,  exit  flow  angle, 
choking,  critical  and  supercritical  flow  conditions.  In  this 
paper,  a  further  improvement  of  the  method  is  attempted. 

1.  Assumptions,  basic  equations,  and  boundary  conditions. 

Assumptions:  perfect  gas,  absolutely  steady,  or  rela¬ 

tively  steady  Insulated  from  the  surrounding. 

Basic  equations  (in  relative  coordinate  system): 

Continuity  equation  v • (pWO  —  0 

Momentum  equation  —'Vp/p-*dW/dt—u>lr+2otxW 

Energy  equation  rf/ -•  0 

Equation  of  state  p  —  pRT 

Thermodynamic  relation  Tdt  —  dh  —  dp/P 


Expression  for  loss  As  —  /fin  {1  +  (y,  4-  Y,  +  y^Ml  —  (1  +  (*  —  1)A/J/2)‘^TH 

I 

Expression  for  exit  flew  ancles  fi,  —  at  —  ^  c,A/r'  +  *t/c 

i  - 1 

This  sot  of  simultaneous  equations  can  be  used  to  solve 
the  flow  parameters  at  the  subcritical  conditions.  The  critical 
and  supercritical  cases  will  be  discussed  later. 

The  boundary  conditions  include  the  dimensions  of  the  flow 
channel,  the  blade  ancle  at  the  exit  of  the  blade,  the  stagna¬ 
tion  pressure  and  stagnation  temperature  of  the  flow  at  the 
inlet,  the  flow  rate  at  the  inlet  or  the  back  pressure  at  the 
exit.  To  solve  for  the  losses,  some  characteristic  dimensions 
of  the  blade  must  also  be  given. 


Figure  1.  The  enthalpy-entropy  relation  of  the  turbine  stage. 
2.  The  losses  and  the  exit  flow  angle 

The  coefficients  of  pressure  loss  are 

Yh  -  (p?  -  p?)/(p?  -  />,)  y,  -  ( pZ '  -  pTVipr  -  p>\ 

The  relations  among  the  parameter's  are  shown  in  Figure  1. 


i 
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The  coefficients  of  stagnation  pressure  loss  can  be 
expressed  by  the  following  equations: 

<tn  «=  1/{1  4-  yw[l  -  (1  +  (*  -  l)A/i/2)T^n 
«■«  —  1/{(I  +  Y.[l  -  (1  -1-  (k  -  1)A/J/2)^1} 

The  stagnation  efficiency  of  the  turbine 

Vt  -  (//j  -  “  «i) 

The  subcritical  and  supercritical  cases  should  be  consid¬ 
ered  separately  for  the  calculation  of  losses.  For  the  sub- 
critical  condition,  the  losses  of  the  flow  past  the  blades  can 
be  classified  into  the  profile  loss,  secondary  flow  loss,  and 
leakage  loss.  The  following  expressions  can  be  applied  for 
common  blade  profiles: 

Profile  loss 

y,  •=  —  -4,  if 'I  ty,<.^»  +  («i/&)!(Yf<#,*’«i>  —  Vf(«l«tn)K5CB1,//>) 

/  "  I  ' 

(  when  '*  <  —  3° ) 

v*  -  (1  +  i  “••'-)  I  '<>*  o,  +  («,/&)*( yft,,-w  -  yK-l-»)K5c«1I/^)-.». 

t  M  7 

(  when  *«  ^  —  3°) 


Secondary  flow  loss 

v,  —  2 l(s,:c)b  cosVj(  <g  £1  —  >g0j)'  *  \A  +  Og0t  +  'g ft)V(A  cos/?,) 

Dunham  and  Came  suggested  that  «*  0.0334 

can  be  adopted  to  account  for  the  boundary  layer  effect  at  the 
end  walls. 

» 

Leakage  loss  y4  —  2H6k  co»'/f;(  ig£,  —  *g/?i)J  *  v^4  +  ( ‘g/^  +  h 


blade 


The  coefficient  of  stagnation  pressure  loss  of  the  guide 

?  **-1/{l  +  (V,  +  V,  +  y.)ll  -u  +(*-l)Mi/2)'r?r1) 
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The  efficiency  of  the  rotor  blade 

i,  -  (//2  -  lit)/ (iii  -  //,n  +  ( y r  +  y,  +  y.)  • 

11  -  (i  +  (k  -  DAU./zp^ir^) 

Entropy  increase  of  the  guiding  and  rotor  blades: 

Ajn=“-  KIntfv  A/k  =  Ik/(1  —  *•  > ]  •  AMn{[//2  +  qT(l'ttU}  —  Vg,«t)/J]/ 

i  nt  h-  ov<.  -  f»«2')/;i) 

To  obtain  the  parameters  for  which  the  losses  are  taken 
into  account,  the  following  set  of  simultaneous  equations 
must  be  solved. 

Continuity  equation  upstream  of  the  throat 

i»« 

C  ■=■  U  +  (k  -  I ) A/ i/2 1  J"-*>  M.  •  C.f"*'- 
Continuity  equation  downstream  of  the  throat 
F./CF, cosft)  —  {[1  +  (k  -  l)Mi/2]/[l  +  (s  -  l)A/]/2]}'5^-(M2/Af.)e-A'"* 

Critical  condition  of  the  flow  rate  dG/dM.^0 

Momentum  equation 

+  (*  _  \  )M\ny&(FjF.) 

CO*  P.  -  11  +  (*  -  l)A/l/2]"i^“"  +  (stfcA/.)[l  +  (k  -  l)Afi/2]1^ 

•  [M.  -  M,y/l  +  (*  -  1  )Af  i/2  cos  (ft  -  ft)A/l  +  (*  -  l)Mj/2  ] 

Expression  for  less 

a*.  -  (/»//)/?  in{  i  +(y,  +  y<  +  y«)  • 

ll .  —  (l  ■+(«  — i)a#!/2)t^1) 

For  the  supercritical  condition,  another  set  of  simultan 
eous  equations  must  be  used.  For  the  triangular  region  at 
the  tail  (see  Figure  2),  the  following  equations  can  be  given 

Continuity  equation 

F«/(Fj  cos  ft)  -  (U  +  (*  -  O  A/ 1/21/ 

[1  +  (*  -  l)Afj/2l}'^(A/J/A/.)f-A'« 
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Figure  2.  Triangular  region  at  the  tail  of  the  turbine 
cascade . 

Momentum  eauation 

+  (*  -  1)A/J/21^(F,/F.)  cos  ft 

-  [1  +  (*  - 

+  *^A/.tl  +  (*  -  l)A/;/2p^[M.  -  Aft  • 

Thermodynamic  relation 

•  Vl  +  (*  -  1)  A/i/2  cos  (ft  -  ft)/\/l  +(*-l)A#J/2  ] 
e"A'"»  —  P°/p° 

For  the  supercritical  case,  this  should  be  satisfied: 

A»,  «-  A/.i 

Combined  with  the  simultaneous  equations  derived  earlier, 
(except  the  subcritical  losses  and  the  expression  for  the 
suberitieal  exit  flow  angle),  these  equations  can  be  solved 
to  obtain  the  flow  field  solution  under  supercritical  condi¬ 
tions.  During  the  solution  process,  the  nonsynchronism  of 
different  parameters  may  lead  to  false  choking.  Hence  the 
synchronized  iterative  method  is  used  in  this  paper  to  avoid 
the  occurrence  of  false  choking.  In  the  program,  the  simul¬ 
taneous  equations  are  first  discretized,  and  then  converted 
into  the  finite  difference  forms  at  the  defined  mesh  points. 
The  relaxed  iteration  method  is  employed  to  solve  this  set  of 
equations  and  the  flow  parameters  at  different  intervals  are 
obtained,  ensuring  the  synchronism  of  different  parameters 
during  the  computation.  The  newest  streamline  pattern  and  the 
corresponding  enthalpy  and  entropy  are  further  obtained,  and 
the  simultaneous  equations  are  solved  again.  This  iteration 
process  is  repeated  until  a  certain  degree  of  accuracy  of  the 
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streamline  position  is  sat.i  sfied  .  Since  the  eoef  f  i  c ;  enis  of 
lor. see.  and  exit  flow  angle  do  not  have  to  be  given  for  each 
deal r.n  condition  In  ♦'his  calculation  process,  it  Is  possible 
to  conduct  the  c amputations  for  different  off-desir.n  condi¬ 
tions  continuously.  Hence  the  computing  efficiency  is  greatly 
increa  red . 

III.  CALCULATED  RESULTS  AMD  ANALYSIS 

Based  on  the  set  of  simultaneous  equations,  a  program  is 
developed  and  Is  applied  to  the  calculation  of  tens  of  exam¬ 
ples  and  design  cases.  Fairly  satisfactory  results  are  obtained. 
The  convergence  speed  of  the  calculation  is  fast.  Only  a  few 
iterations  of  the  streamlines  are  sufficient  to  reach  a  1/1000 
accuracy  for  the  relative  location  of  the  streamlines.  The 
results  obtained  from  the  supercritical  calculation  are  also 
reasonable,  and  are  found  to  be  in  agreement  with  the  experi¬ 
mental  results. 

Calculations  are  conducted  for  the  turbine  of  the  250  kw 
gas  engine  developed  partly  by  the  Institute  of  Mechanics,  and 
the  turbine  of  the  foreign  JT3D-3B  engine.  The  results  are 
satisfactory,  as  described  below: 

1.  Calculated  results  for  the  250kw  turbine 

This  is  a  two-st.age  turbine,  designed  according  to  con¬ 
stant  operating  conditions  along  the  blade  span.  The  blade 
profile  is  one  of  the  commonly  used  airfoil  profiles.  Single- 
stage  gas  turbine  tests  have  been  conducted  for  the  first 
stage,  and  the  second  stage  separately.  Experimental  data  on 
the  second  stage  are  more  complete.  Hence  the  experimental 
results,  and  the  calculated  results  of  the  second  stage  are 
compared  here. 
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shown  In  Figure 


The  flow  channel  in  live  i’hOkw  tnrbJ  ne  in 
3.  Figure  <1  shown  the  relation  of  the  expansion  ratio 
and  the  reduced  flow  rate  c,<Jr° /p°  for  th°  entire  off-design 
ran,te  of  the  test.  The  solid  line  in  the  figure  indicates 
the  calculated  result  obtained  b y  the  present  method,  while 
the  circles  indicate  the  experimental  results.  It  is  apparent 
that  the  calculated  results  agree  well  with  the  experimental 
results.  The  error  is  within  1?. 


Figure  3.  Flow  channel  in  the  250kw  turbine. 

1  -  air  flow 

From  Figure  *4,  it  can  also  be  observed  that  close  to  the 
choking  condition,  both  the  experimental  data  and  the  theoreti¬ 
cal  results  from  the  calculation  approach  a  horizontal  line 
(that  is,  the  limit  of  the  choking  flow  rate)  with  good  con¬ 
sistency.  The  relation  between  the  turbine  efficiency  nT 
and  the  expansion  ratio  tt^  is  shown  in  Figure  5.  It  is  obvious 
that  the  calculated  results  are  close  to  the  experimental 
results  for  the  entire  test  range.  The  error  ’is  only  about 
155.  Figures  6  and  7  compare  the  distribution  of  the  flow  para¬ 
meters  in  the  radial  direction.  The  radial  distribution 
patterns  of  the  flow  angle  and  the  Mach  number  at  the  exit  plane 
of  the  guiding  blade  are  shown.  There  is  a  significant  differ- 
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orcnro  i>f «^n  l  h<>  r:>.  1  <■•••}  ->nd  experiment  n  1  results  at 

certain  indi  vidua  1  noin' s  only.  They  are  fairly  consistent 
in  the  main si  ream  region.  The  calculated  results  for  the 
JT3D-3B  are  presented  ir.  Table  1. 


Figure  4.  The  relation  between  and  tt-,  of  the  250kw 

II  stage.  e 


Figure  5.  The  relation  between  hm  and  it-  of  the  250kw-II 


I 
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Figure  6.  Radial  distribution  Figure  7.  Radial  distri¬ 
bution  of  bution  of 


Table  1.  Comparison  of  the  parameters  and  calculated  results 
of  the  JT3D-3B  high  pressure  turbine 


Item 

Given  Parameters 

.Calculated  Result 

stagnation  pressure  at 
inlet  P° 

13.182  kg/cm3 

13.182  kg/cm3 

stagnation  temperature 
at  inlet  T° 

1160.6  K 

1160.6  K 

intake  mass  flow  rate  G 

88.335  kg/sec 

88.01  kg/sec 

rotational  speed  n 

9800  rpm 

9800  rpm 

stagnation  pressure  at 
exit  P^ 

6.151  kg/cm2 

6.086  kg/cm 2 

stagnation  temperature 
at  exit  T^ 

984.9  K 

982.7  K 

expansion  ratio 

2.143 

2.166 

stagnation  temperature 
ratio  6 

1.178 

1.181 

We  can  now  conclude  the  following: 

(1)  This  method  is  suitable  for  calculating  the  direct 
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problem  of  transonic  turbines  at  both  the  design  point,  and  the 
off-design  condition.  The  losses  and  the  exit  flow  ancles 
can  be  obtained.  The  choking  flow  condition  can  also  be  cal¬ 
culated.  There  is  no  difficulty  in  solving  the  flow  parameters 
for  the  critical  and  supercritical  cases. 

(2)  The  losses  of  common  blade  profiles  (e.g.  the  blade 
profile  of  the  250  kw  and  JT3D-3B  turbines)  calculated  by  this 
method  are  comparatively  reasonable.  A  comparison  between  the 
calculated  and  the  experimental  result  for  the  entire  test 
range  shows  a  very  satisfactory  consistency. 

(3)  The  agreement  between  the  experimental  results  arid 
the  calculated  variation  of  reduced  flow  rate  versus  expansion 
ratio  close  to  the  choking  conditions  indicates  that  relatively 
reasonable  results  can  be  obtained  for  the  throat  choking  flow 
rate  by  using  the  critical  and  supercritical  conditions  with 
losses  taken  into  account. 

(4)  The  simultaneous  equations  including  the  expression 
for  losses  and  exit  flow  angle  are  applied  in  this  paper  to 
solve  for  the  flow  parameters  of  the  turbine.  Hence  the  flow 
parameters,  losses,  and  exit  flow  angles  can  be  obtained  at 
the  same  time.  This  approach  enables  us  to  find  a  reasonable 
choking  condition  quite  easily.  It  also  makes  the  continuous 
computation  for  off--design  conditions  possible. 

(5)  The  consistency  of  the  experimental  and  the  calcula¬ 
ted  results  of  the  relation  between  the  reduced  flow  rate  and 
the  expansion  ratio  suggests  that  the  computation  method  is 
appropriate  for  the  exit  flow  angle  of  common  blade  profiles 
within  the  test  range.  This  conclusion  can  be  observed  from 
Figure  6.  The  exit  flow  angle  along  the  blade  height  is  also 
consistent . 


(6)  Certain  ar.cun.pt  1  on:;  are  made  iri  this  method.  That 
is  why  the  limiting  Mach  number  is  restricted  by  the  limiting 
design  conditions  of  the  turbine.  The  error  will  increase 
if  the  Mach  number  overshoots  the  limit.  Simplifications  are 
also  made  for  the  boundary  layer  loss  at  the  end  walls. 

Hence  there  is  a  significant  error  of  the  flow  parameters  near 
the  end  walls.  This  phenomenon  requires  further  investigation. 
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empirical  coefficient 
empirical  coefficient 
chord  length 

constant,  or  blade  thickness 
radius  of  curvature  of  the 
tail  of  the  blade  back 
area 

flow  rate 

absolute  stacnation  enthalpy 
blade  height  or  static  enthalpy 
relative  stagnation  enthalpy 
length  of  the  central  stream¬ 
line  of  the  channel  between 
blades 

temperature 

time  or  distance  between 
blades 

circumferential  velocity 

absolute  velocity 

relative  velocity 

loss  coefficient 

profile  angle  of  the  guide 

blade 

profile  angle  of  the  rotor 
blade  or  the  flow  angle  of 
the  guide  blade 


u>  angular  velocity 

gap  at,  the  tip  of  the  blade 
Subscri pt s : 

1  guide  blade  input 

2  guide  blade  exit 

3  rotor  blade  exit 
k  leakage 

N  guide  blade 

n  throat  or  reduced  quantity 
nt  from  the  throat  to  the  exit 
ne  critical 
0  tangential 

p  blade  profile 
r  radial 
R  rotor  blade' 
s  secondary 
T  turbine 
W  relative 
z  axial 

Superscripts: 
o  stagnation  state 
o'  relative  stagnation  state 


density 

angle  of  attack 

recovery  coefficient  of  stagnation  pressure 


efficiency 
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Abstract 

Wc  obt«in  a  group  of  simultaneous  equations  including  equations  of  S,  stream  sur¬ 
face,  simplified  S ,  stream  surface  (for  use  in  the  solution  of  the  triangle  region  at  tail 
of  cascade),  flow  loss  and  deviation  angle  correlations  for  calculation  of  direct  problem 
in  turbomachino.  By  the  solution  of  these  equations,  the  design  point  as  well  as  off- 
design  point  performances  can  be  calculated. 

The  essential  difference  in  treating  sub-critical,  critical  and  super  critical  flow 
condition  is  considered.  Influence  of  losses  is  taken  into  account  in  the  calculation  of 
choking  conditions.  We  adopt  a  unified  approach  to  solve  these  equations.  Losses, 
deviation  angles  and  other  flow  parameters  are  obtained  in  the  same  iterative  process, 
so  that  false  choking  condition  will  not  occur  and  convergence  will  be  rapid.  Computer 
programm  was  then  completed  and  may  be  used  in  the  calculation  of  the  off-design 
performance  of  transonic  turbines. 

In  applying  this  programm  to  the  calculation  of  the  performance  of  exist  turbine 
stage,  we  have  got  satisfactory  results.  The  reduced  flow  rate,  efficiency  and  expan¬ 
sion  ratio  are  fairly  consistent  with  the  experimental  results.  The  difference  between 
the  computing  results  and  experimental  data  is  within 
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